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ABSTRACT 


A  computer  simulation  of  the  coverage  of  satellite  borne  radars 
with  mechanical  and/or  electronic  scanning  mode  has  been  implemented  in  APL 
language.  The  number  of  revisits  and  the  revisit  time  are  calculated  for  a 
set  of  points  of  reference.  The  programs  are  producing  either  paper  print¬ 
out  or  graphic  displays  on  a  Tektronic  screen.  This  simulation  is  a 
component  of  a  software  package  designed  to  stucjy  space  based  radar 
surveillance  systems. 


RESUME 


On  presente  une  simulation  par  ordinateur  en  language  APL  de  la 
couverture  de  radars  montes  sur  satellites  a  balayage  electronique  et/ou 
mecanique.  Les  programes  produisent  une  sortie  imprimee  ou  afflichent  les 
donnees  sur  ecran  Tektronix.  On  calcule  1'intervalle  entre  les  visites  et 
le  nombre  de  revi sites  pour  un  ensemble  de  points  de  reference.  Cette 
simulation  fait  partie  d'un  ensemble  de  logiciels  destines  a  l'analyse  de 
systemes  de  surveillance  par  radars  montes  sur  satellites. 
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Introduction 

TSATCO  and  SATCO  are  two  complementary  workspaces  whose  function 
is  to  calculate,  for  a  set  of  reference  points,  the  number  of  revisit  and 
the  revisit  time  associated  to  a  particular  system  of  satellites  described 
by  the  variable  SYSAT.  TSATCO  generates  a  graphic  display  on  a  Tektronix 
screen  and  SATCO  generates  a  print-out  of  a  few  variables  of  interest. 

TSATCO  is  used  to  check  the  results  obtained  from  SATCO  and  to  set-up  the 
time  interval  that  is  appropriate  to  a  particular  situation.  The  components 
of  the  system  of  satellites  are  defined  by  the  user  through  the  definition 
of  the  variable  SYSAT. 

Function  of  the  Workspace  TSATCO  and  SATCO: 

-  Calculates  the  position  of  the  satellites, 

-  Calculates  the  coverage  area  of  the  satellites, 

-  Checks  if  one  or  more  of  the  check  points  are  within  the  coverage 
areas, 

-  Computes  the  revisit  times  and  the  number  of  revisits,  for  every 
reference  point, 


-  TSATCO  displays  the  satellite  position  and  the  coverage  area  if  the 
coverage  area  includes  one  of  the  reference  points, 

-  SATCO  prints  a  few  parameters  of  interest, 

-  SATCO,  at  the  end  of  one  execution  of  CONTROL,  is  saved  in  a 
workspace  named  DATACO.  The  values  of  the  revisit  time  and  of  the 
number  of  revisits  are  thus  saved  and  it  is  possible  to  further 
process  them,  before  another  execution  of  CONTROL,  with  the  program 
COURBE  and  HISTO. 


Mode  of  Operation 


The  calculations  are  made  at  a  regular  time  interval  from  a  time 
origin,  both  are  set  by  the  operator.  The  calculations  are  separated  in 
pockets  of  a  certain  number  of  time  points  such  that  the  parallel  processing 
capabilities  of  APL  are  used  to  the  maximum  without  exceeding  the  available 
memory  of  the  workspace.  The  number  of  time  points  processed  in  parallel  is 
given  by  the  variable  NLIM  that  the  user  can  change  at  will  on  line  22  and 
21  of  TSATCO  and  SATCO  respectively. 


Reference  points: 


The  reference  points  are  located  at  longitude  0  and  latitude  90°, 
80°,  70°,  60°,  50°,  40°,  30°,  20°,  10°  and  0°  in  geographical  coordinates 
(EGR  system). 


Example  of  Results 

We  included  an  example  of  the  output  of  TSATCO  (see  Fig.  1)  and 
SATCO  (see  Table  1)  with 

-  starting  time  of  the  simulation:  0 

-  time  increment:  5 

-  duration  of  the  simulation:  150 

We  had  a  system  of  two  satellites  on  the  same  circular  polar  orbit 
at  an  altitude  of  300  km  but  spaced  in  time  by  15  minutes.  One  of  the 
satellites  was  of  the  mechanically  scanning  type  and  the  other  was  of  the 
electronically  scanning  type.  The  SYSAT  array  was: 
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List  of  the  systems  of  coordinates  used  in  the  programs 


Name 

Center 

Type 

Rotation^3  ^ 
with  the 
Earth 

Orientation  of  the  axis 

ERF 

Earth 

rectangular 

No 

Z-North  Pole;  X  at  0°  lat. 
and  0°  long,  at  t=0;  see 
Fig.  2. 

ESF 

Earth 

spherical  Z^ 

No 

from  ERF ;  see  Fig.  2. 

ERR 

Earth 

rectangular 

Yes 

Z-North  Pole;  X  at  0°  lat. 
and  0°  long;  at  anytime 

ESR 

Earth 

spherical 

Yes 

From  ERR,  see  Fig.  2. 

ERE 

Earth 

rectangular 

No 

Z-perpendicular  to  the 
orbital  plane;  X  toward 
perigee;  see  Fig.  3 

ESE 

Earth 

spherical  Z^ 

No 

From  ERE;  see  Fig.  2  and  3. 

EGR 

Earth 

geographic 

Yes 

Altitude:  from  the  center 
of  the  Earth;  longitude, 
latitude;  see  Fig.  4 

SRT 

satellite 

rectangular 

- 

topocentric  (4);  (see 

Fig.  3). 

SST 

sate  1 1 ite 

spherical  Z^) 

- 

From  SRT 

SXT 

satel 1 i te 

spherical  X^ 

- 

From  SRT  (see  Fig.  5) 

SRE 

satel 1 ite 

rectangular 

Z-perpendicular  to  the 
orbital  plane;  X  in  the 
direction  of  the  movement 
(see  Fig.  5). 

SXE 

satel 1 i te 

spherical  X^ 

- 

From  SRE  (see  Fig.  5). 

ORR 

point  of 
observation 

rectangular 

Yes 

Topocentric  (see  Fig.  6). 

(a)  Yes  means  that  the  system  of  coordinates  rotates  with  the  Earth.  No 
means  that  the  system  of  coordinates  does  not  rotate.  It  is  just  tied 
to  the  center  of  the  Earth  and  assumed  to  keep  a  constant  orientation 
relative  to  the  stars.  The  movement  of  the  Earth  around  the  Sun  are 
thus  neglected. 

(b)  Spherical  coordinates  with  the  angle  <t>  measured  from  the  Z  axis  (see 
Fig.  2). 

(c)  Spherical  coordinates  with  the  angle  <|>  measured  from  the  X  axis  (see 
Fig.  5). 


Figure  3:  System  of  coordinates  used  to  describe 
the  motion  of  the  satellites 


LONGITUDE  0  AT 
GREENWHICH 


LONGITUDE  POSITIVE 
TOWARDS  EAST 


Figure  4:  Geographic  coordinates 


POINT  OF  OBSERVATION 
LONGITUDE  -90° 


LATITUDE  45° 


Figure  6:  Topocentric  system  of  coordinates 


Description  of  the  systems  of  coordinates 


fixed  system  of  rectangular  coordinates  centered  on  the  Earth.  The 
origin  of  the  system  is  tied  to  the  center  of  the  Earth  and  the 
system  of  axis  does  not  rotate  with  the  Earth.  The  Z  axis  is 
pointing  to  the  North  Pole  and  the  X  axis,  at  time  t=0,  is  pointing 
towards  the  point  with  longitude  0  and  latitude  0.  The  movement  of 
the  Earth  around  the  Sun  is  neglected  because  the  simulation  is  to  be 
used  for  short  periods  of  time  going  from  a  few  minutes  to  one  day 
(see  Fig.  2). 

fixed  system  of  spherical  coordinates  where  R,  o  and  4>  are  defined  in 
the  usual  way  from  ERF,  41  is  measured  from  the  Z  axis  (see  Fig.  2). 

rotating  system  of  rectangular  coordinates.  The  origin  of  the  system 
is  tied  to  the  center  of  the  Earth  and  the  system  of  axis  does  rotate 
with  the  Earth.  The  Z  axis  is  pointing  towards  the  point  with 
longitude  0  and  latitude  0  all  the  time. 

rotating  system  of  spherical  coordinates  where  R,  0  and  f  are  defined 
in  the  usual  way  from  ERR,  41  is  measured  from  the  Z  axis  (see 
Fig.  2). 

fixed  system  of  rectangular  coordinates.  The  origin  is  at  the  center 
of  the  Earth,  the  Z  axis  is  perpendicular  to  the  orbital  plane  of  the 
satellite  and  the  X  axis  is  pointing  towards  the  perigee.  The 
satellite  turns  in  the  positive  direction. 

fixed  system  of  spherical  coordinates.  Where  R,  0  and  4>  are  defined 
in  the  usual  way  from  ERE.  4>  is  measured  from  the  Z  axis  and  0  is 
the  true  anomaly  (see  Figs.  3  and  10). 

rotating  system  of  geographic  coordinates.  The  origin  is  tied  to  the 
center  of  the  Earth,  the  radial  distance  is  measured  from  the  origin. 
The  latitude  is  positive  in  the  Northern  Hemisphere  and  negative  in 
the  Southern  Hemisphere.  The  positive  direction  for  the  longitude  is 
Eastward  (see  Fig.  4). 

rectangular  and  topocentric  system  of  coordinates  centered  on  the 
satellite.  The  X  axis  is  pointing  away  from  the  center  of  the  Earth 
and  the  Z  axis  is  pointing  Northward  (see  Fig.  3  and  6). 

spherical  system  of  coordinates  centered  on  the  satellite  where  R,  0 
and  <4  are  defined  in  the  usual  way  from  SRT.  4>  is  measured  from  the 
Z  axis  (see  Fig.  1). 

spherical  system  of  coordinates  centered  on  the  satellite  where  R,  0 
and  4>  are  defined  from  SRT,  according  to  Fig.  5.  41  is  measured  from 

the  X  axis. 

rectangular  system  of  coordinates  centered  on  the  satellite.  Z  is 
perpendicular  to  the  orbital  plane  and  Y  is  tangent  to  the  orbit  (see 
Fig.  7). 


spherical  system  of  coordinates  centered  on  the  satellite  where  R,  j 
and  <p  are  defined  from  SRE  according  to  Fig.  5.  <f>  is  measured  from 
the  X  axis. 


rotating  topocentric  rectangular  system  of  coordinates.  The  origin 
of  the  system  is  at  the  point  of  observation.  The  X  axis  is  pointing 
away  from  the  center  of  the  Earth  and  the  Z  axis  is  pointing 
Northward  (see  Fig.  6). 


The  units  used  in  that  workspace  are: 

-  for  the  time:  minutes 

-  for  distance:  kilometers 

-  for  angle:  degree 

-  for  speed:  ki lometers/mi nute. 


Description  of  the  input  variables 


contains  the  parameters  of  the  satellite  orbits  and  the  type  of 
radar  set  carried  by  the  satellite.  The  format  of  SYSAT  is 
illustrated  in  Figure  8.  The  type  of  scanning  performed  by  the 
radar  is  set  in  column  10:  0  stands  for  electronic  scanning  and  1 

stands  for  mechanical  scanning  (see  Fig.  8). 
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Documentation  of  the  programs 


The  documentation  of  the  programs  contains,  for  each  program,  the 
list  of  the  input  variable,  their  dimension,  their  nature,  their  source 
program,  system  of  coordinates  and  units.  A  list  of  the  output  variables, 
their  dimensions,  their  nature,  their  destination  program,  the  system  of 
coordinates  and  their  units  is  also  provided.  It  is  also  stated,  for  every 
variable,  between  the  parenthesis  following  its  name,  if  the  variable  is 
global  (G),  local  dumny  (LD)  or  local  result  (LR). 


f  unction:  calculates  the  polar  angle  associated  to  a  vector  who  has 
certain  component  X  and  Y. 


Var I  ah les 


Output  Dimensions 
Varl ables 


Nature 


Destination  System  of  Units 
program  Coordl  nates 


number 

polar  angle  associated  to  the 

ANGROT 

components  X  and  Y  of  a  vector 

TVRCSO 

-  21 


m 


K  ■' 


V  J-  »'  1 1I  '.I  n  V  w '  .  II  U  ■■  U/  Uv-  Tw  I"  I  l  »■  i  in  !-»■—{- 


ANGRV 


Function:  calculates  the  polar  angle  associated  1o  each 
vector  having  certain  components  X  and  Y. 


1  nput 
Varl ables 

Dimensions 

Nature 

Source 

program 

System  of 
Coordl  nates 

Units 

AX(LO) 

vector 

contains  the  X  corrponents  of  the 
vectors 

TRECSPX 

rectangular 

unitless 

AY(LO) 

vector 

contains  the  Y  components  of  tne 
vectors 

TRECSPX 

rectangu  1  ar 

un 1 1 1  ess 

1 


1 


i 


ANGROT 


Function: 


calculates  the  three  rotations  around  the  axis  Y,  X  and  Z  that 
are  needed  to  go  from  a  system  XYZ  to  a  system  X'Y'Z' 

(see  Fig.  9).  The  rotation  matrix  necessary  to  make  the  trans¬ 
formation  and  the  Inusrse  transformation  are  calculated. 


1  nput 
Variables 

D 1  mens  Ions 

Nature 

Source 

program 

System  of 
Coordl  nates 

Units  r  f- 

VZ(LD) 

vector  ( 3 ) 

coordinates  of  Z'  In  the  XYZ 
system 

* 

rectangular 

unitless 

i 

VX(LO) 

vector  (3) 

coordl  nates  of  X'  In  the  XYZ 
system 

* 

rectangu  lar 

unit  less 

-  -1 

AREP(G) 

number 

angle  associated  to  a  vector 
having  a  certain  X  and  Y 
coordinate. 

ANGR 

— 

radians 

11 

RX  (G ) 

3x3 

Matrix  of  rotation  that  executes 

a  rotation  of  PHX  around  the  X 

axl  s. 

RAX 

rectangular 

unit  less 

"  i 

RY  (G) 

3x3 

Matrix  of  rotation  that  executes 

a  rotation  of  PHY  around  the  Y 

axl  s. 

RAY 

rectangular 

unitless  ~ 

• 

3x3 


Matrix  of  rotation  that  executes 
a  rotation  of  PHZ  around  the  Z 
axl  s. 


RAZ 


rectangular 


unit  less 


ANGROT 


Outpu  t 

Var 1  ah les 

Dimensions 

Nature 

Destination 

program 

System  of 
Coordl  nates 

Units 

MR  (G ) 

(3  x  3) 

matrix  of  rotation  that  trans¬ 
forms  X'Y'Z'  Into  XYZ 

« 

rectangular 

uni t l« 

IMR(G) 

(3  x  3) 

matrix  of  rotation  that  trans¬ 
forms  XYZ  Into  X'Y'Z* 

* 

rectangu  1  ar 

un  1 1  le 

BROWNS 

Function:  calculates  the  angle  f  (the  true  anomoly,  see  Fig.  10) 

associated  to  an  Input  vector  containing  various  time  for 
a  cental n  orbit. 


Variables 


Dimensions 

Nature 

Source 

program 

vector  (NT ) 

contains  the  time 

POSSAT 

vector(NS ) 

eccentricity  of  the  orbit 

POSSAT 

number 

number  of  time  point 

POSSAT 

number 

number  of  satellite 

POSSAT 

BROWNS 

Output  Dimensions 

Nature 

Destination 

System  of 

Variables 

program 

Coord 1  nates 

m 

y  f*(LR) 

vector  (NT ) 

true  anomoly 

POSSAT 

Is  given  by  a  trascendenta  1  equation:  the  Brown's  method, 
as  described  In  reference  1,  has  been  used  to  solve  T. 


COURSE 


Function:  draws  the  curve  specified  by  the  vector  VV. 


Function: 

draws  a  circle  that  represents  the  Earth. 
Input  and  no  output  variables* 

That  program  has  no 

1  nput 

01  mens  Ions  Nature 

Source  System  of 

Variables 

program  Coordinates 

h 

(10  v  3) 

contains  the  location  of  the 
reference  points 

CONTROL 

ESR 

O] 


1 


" jnct 1  on: 

calculates 

the  matrix  of  rotatl 

on  that 

transforms 

a  ERF  I nto  a  SRT  or 

vl  ce  versa 

1  npjt 

D I tensions 

Nature 

I  Source 

System  of 

V ar 1 ab les 

program 

Coordl  nates 

PO(LD) 

vector  (3) 

position  of  the  origin  of  the 
SRT  system  In  the  ERF  system 

CONTROL 

POSSAT 

ERF 

km. 

IMR(G) 

(3  x  3) 

matrix  of  rotation  to  transform 

ERF  Into  SRT 

ANGROT 

unitless 

D I  mens  Ions 


Variables 


Nature 


Destination  System  of 
program  Coordinates 


MRO(LR>  I  (3x3)  matrix  of  rotation  to  transform  CONTROL 
ERF  I nto  SRT  or  ERF  I nto  ORR 


unit  less 


(3  x  3) 


matrix  of  rotation  to  transform 
SRT  I nto  ERF 


POSSAT 


unitless 


POSSAT 


Function:  calculates  the  position  of  the  sate  I  I  Ite  I  n  the  ERF  and  the  ORR 
system  of  coordinate.  It  also  calculates  translations  and 
matrix  of  rotation  to  be  used  In  later  calculations. 


I  nput 
Varl ables 


PT(LD) 


SYSAT(G) 


F«5> 

MR  (G ) 


MR  (G ) 


MR(G) 


MRO(G) 


XDR.YDR 

ZDR(G) 


D  I  mens  Ions 


vector!  3) 


(NS,  11) 


(NS.NT) 
(3  x  3) 


(3  x  3) 


(3  x  3) 


(NT  x3x3) 


number 


Nature 


Source 

program 


vector  containing  the  starting 
time,  the  time  Interval  and  the 
termination  time  of  the  slmula- 
t  Ion 


CONTROL 


contains  the  parameters  of  the 
orbit  of  the  satellites 


user 
def I ned 


true  anomoly 


matrix  to  go  from  the  ERE  to  the| 
ERF  system  (line  12) 


BROVINS 

ANGROT 


matrix  to  go  from  the  SXT  to  the| 
ERF  system  (line  13) 


0P0 


matrix  to  go  from  the  SXT  to  the| 
ERE  system  (line  14) 


ANGROT 


matrix  to  go  from  the  ERF  to  thej 
ORR  system 


TOPO 


translation  to  go  from  ERF  to 
ORR 


TOPO 


System  of 
Coord!  nates 


POSSAT 


# 


Outpu  t 

Variables 

P  1  mens  1  ons 

Nature 

Destination 

program 

System  of 
Coordl  nates 

MRO(G) 

(NT, 3,3) 

matrix  to  go  from  the  ERP  to 
the  ORR  system. 

POSSAT 

TDE,  TOM 

- 

MTR(G) 

(NT, NS,  3) 

translation  to  go  from  ERE  to 
ORR 

POSSAT 

TDE,  TDM, 
TFM,  TFE 

PSG(G) 

(3, NS, NT) 

satel llte  position  1  n  the  ERE 
system 

CONTROL 

- 

PSO(G) 

(3, NS, NT ) 

satel 1 Ite  position  1 n  the  ORR 
system 

CONTROL 

- 

MRS  (G ) 

(NS, NT,  3, 3) 

! 

matrix  to  go  from  the  SXT  to 
the  ERE  system 

TDE,  TDM 

- 

Units 


un i t less 


un 1 1 1  ess 


km 


km 


un 1 1 1  ess 


MRS*:  (!> ) 


(NS,NT,3,3)  matrix  to  go  from  the  SXE  +o 
the  ERE  system 


TDM 


un 1 tless 


I 


s.j  h  i  i, 'urn  r -e  coordinate  of  the  reference 
point-,  >ro m  Mie  rSR  system  to  the  SSX  system 


:  n  i  j  * 

:  1-i.sMes 

I  th<is  1  ons 

Nature 

Source 

program 

System  of 

Co.x'dl  nates 

Units 

/  '  ■) '  1 

prlqln  of  the  3GX  system 
i  (position  of  the  satellite) 

'JON  TP  01. 

■ 

FRF 

km. 

PRS  C.  n ) 

10  x  x 

1 

n  ">s  i  f  i  on  'if 

•Ol  'lf'5 

rw^r 

^  S' 

TIG) 

number 

t  1  T** 

yiNmi. 

ml  n. 

prs(G) 

10  x  3 

nos  1 1 1  on  of  roforence  points 

TSPHREC 

ERR 

km. 

PRS<3) 

10  x  3 

position  of  reference  points 

V0P0C 

SRT 

km. 

PSX(G) 

10x3 

position  of  reference  points 

TRECSPX 

SXT 

km. 

radl  ans 

s  =  POSITION  OF  THE  SATELLITE,  ORIGIN  OF  THE  SYSTEM  OF  COORDINATE 
R.P.  =  REFERENCE  POINT 


Figure  11:  A  reference  point  as  seen  from  the  satellite 
using  the  SXT  or  the  SXE  system  of  coordinates 


runcMon:  to  transform  the  coordinate  of  the  reference 

points  from  the  ESR  system  to  the  33E  system 


I  nput 
Var 1 ables 


v0"  (1.0  i 

prG(l  >> 

TIG) 

PRS(G> 

PRS(G) 

PSX(G) 


Dimensions 

Nature 

Source 

program 

System  of 
Coord!  nates 

_ _ _ _ 

vector  ( 3 ) 

origin  of  the  SSE  system 

CONTROL 

ERF 

(position  of  the  satellite) 

10  x  3 

position  of  reference  points 

CONTROL 

ESP 

number 

time 

CONTROL 

- 

10  x  3 

position  of  reference  pol nts 

TSPHREC 

ERR 

10  x  3 

position  of  ref erence  pol nts 

VOPOE 

SRE 

10  x  3 

position  of  reference  points 

TRECSPX 

SXE 

radl  ans 


Dimensions 


Variables 


Pf)3  (3  > 


vector  (10) 


Nature 


position  of  reference  points 


angle  between  the  nadir  and 
the  az  I mjth  position  of  the 
reference  pot nt  (see  Pig-  11) 


Destl nation 

System  of 

program 

Coordl  nates 

CONTROL 

SXE 

STATO 

CONTROL 

STATM 

SXE 

radl  ans 
radl ans 


jo.;  M  )i: 


stabllsn,  for  the  electronic  redlining  node.  If  the  reference  points  are 
wltila  or  out  of  the  coverage  are.).  The  revisit  time  and  the  number  of  revisits 
are  calculated  when  needed. 


Variables 


PDS(G) 


VGAN(G) 


RDMV  (  3) 


D 1  mens  Ions 

Nature 

Source 

program 

System  of 
Coord!  nates 

_ ______ _ _ _ _ _ _ 

(10,  3) 

position  of  reference  points 

FRSE 

CONTROL 

SXT 

vector  (10) 

angle  between  the  nadir  and 
the  azlnuth  position  of  the 
reference  pol nts 

CONTROL 

PRSE 

SXT 

vector  (10) 

angle  between  the  nadir  and 
the  outside  border  of  the 

coverage  area 

TFE 

SXT 

vector  (10) 

distance  between  the  satellite 

aid  the  plane  containing  the 
outside  harder  of  the  coverage 

area 

TFE 

SRT 

radl  ans 
radl  ans 

radl  ans 


Variables 


NR!  7 '31 


<lV '  (3) 

Rl'V ?  (G ) 
<  V  5(G) 
REV4(G ) 
REVb(G) 
REV6(G ) 
REV  7(G) 
REV8CG ) 
REV9(G) 
REVO(G) 


Dimensions 

Nature 

vector  (ID) 

vector 

each  co'VO'Hnt  contains  the 
number  of  revisits  associated 
to  the  correspond!  ng  reference 
paint. 

revisit  tone  for  reference 

po  1  nt  1 

"  2 

M 

"  3 

Destination  System  of 
program  Coord!  nates 


met!  isfabllsh,  for  the  mechanical  scanning  mode.  If  the  reference  points  are 

->i  I  *"h  I  n  pr  out  of  the  coverage  ar  m.  The  ravlslt  tluva  and  the  numbar  pf  revisit 
are  calculated  when  needed. 


Input  |  Dimensions 
Variables 


Nature 


Source  System  of  Units 

program  Coordinates 


PSX(G)  (10,  3) 


position  of  reference  points  PRSM 

CONTROL 


km.  , 

radl  ans 


PDS(G)  vector  (10)  angle  between  the  nadir  and 

the  azimuth  position  of  the 
reference  pol nts 


CONTROL 


radl  ans 


VGAN(G)  vector  (10)  angle  between  the  nadir  and  TFM 

the  outside  border  of  the 
coverage  area 

ROMAX(G)  vector  (10)  distance  between  the  satellite  TFM 

and  the  plane  containing  the 
outside  border  of  the  coverage 


radl  an  s 


WSCAN(G)  vector  (10)  width  of  the  scan  line 


defined  Ini 
SYSAT 


degrees 


D I  mens  Ions 


Nature 


Varl ables 


Destination  System  o‘ 
program  CoordinaK 


NREV(G) 


RFV 1(G) 

RFV2(G) 
REV3(G ) 
REV  4(G) 
REVb(G) 
(  V6< G) 
REV  7(G ) 
REV8(G> 
REV9(G ) 
REV0(G) 


vector  (10)  each  component  contains  the 

number  of  revisits  associated 
to  the  corresponding  reference 
pol  nt. 

vector  revisit  time  for  reference 

point  1 


SYSVEC 


''•motion.  I-)  calculate  f'jn  an  Initial  vector  a  system  of  vectors  located 
at  ‘I  Juj.-ij,  from  each  other.  That  routine  inay  be  used  to  find 
the  parameters  to  be  put  In  SY3AT  as  direction  cos  He  of  the  normal 
to  the  orbl t. 


1  npu  t 
Var 1 ab les 

0  1  mansions 

Nature 

Source 

program 

System  of 
Coord!  nates 

7(1. C) 

vncfor  C5) 

Initial  vector 

ERF 

NC.b) 

number 

number  of  vectors 

VIC.'i  i 

n.|.nb»*r 

jujI  }  between  the  vectors 
(degrees) 

SYSVFC 


hitput 
V  v  i  ab les 


S(G) 


Cl  mans  Ions 


Nature 


pestl  nation 
program 


System  of 
Coord!  nates 


N  X  3 


each  row  contains  one  of  the 
new  vectors 


ERE 


■nof  I  on: 


cal  culates  the  pos  I  H  >n  of  Hi-*  *•  .**••.  of  the  coverage  area 
and  displays  them  (electronic  scanning)  See  Figures  '2  and  13. 


1  nput 
Variables 

D I  mens  1 ons 

Nature 

Source 

program 

ir  (0) 

n umber 

distance  from  the  satellite  to 

the  center  of  the  earth 

TFE 

CONTROL 

XRC  <G ) 

vector  ( 2 ) 

distance  between  the  center  of 
the  Tarth  and  the  planes 
containing  the  boundaries  of 
the  coverage  area 

TFE 

CONTROL 

RRC(G) 

vector  ( 2 ) 

radius  of  the  boundaries  of 
the  coverage  area 

TFE 

CONTROL 

MRS(G) 

NSxNTx3x3 

matrix  to  go  from  the  CRT  to 
the  ERF  system 

possat 

MRO(G) 

NT  x3x3 

matrix  to  go  from  the  ERF  to 
the  ORR  system 

POSSAT 

TRSO(G) 

vector  (3) 

translation  to  go  from  a  SRT 
to  a  ORR  system 

TFE 

;  r 

n jmber 

Index  of  the  satellite 

CONTROL 

IN 

number 

Index  of  the  time  point 

CONTROL 

System  of 
Coord! nates 


44 


SATELLITE 

POSITION 


Figure  12:  Parameters  used  in  the  calculation  of  the  coverage  area 
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DM 


rjnctlon:  calculates  the  position  of  the  borders  of  the  coverage  area 

and  displays  them  (mechanical  scanning)  (see  figure  14). 


1  nput 
Var 1 ables 

Dimensions 

II 

Nature 

Source  j 

program  i 

Sy  s  t  o  f  | 

Cov  >  j 

U  n  1  +  s 

■VSCA‘1 

number  j 

angular  width  of  t'ie  scan 

^  "  1 

use'  def. 

in  >  Y  3  A  T 

i 

XRC*  (3 ) 

vector  (2) 

distance  between  the  center  of 

the  earth  and  the  planes 
containing  the  boundaries  of 
the  coverage  area* 

TFM 

5Jr 

<■  T. 

RRC*(G) 

vector  (2) 

radius  of  the  boundaries  of 
the  coverage  area 

TFM 

- 

km. 

MRSE(G) 

NSxNT  x3x3 

matrix  to  go  from  the  SRF  to 
the  FRF  system 

0TV,AT 

l 

■•T(0(  ',) 

NT<3x3 

matrix  to  go  from  the  FRF  to 
the  ORR  system 

* 

dso<:;> 

vect  or  ( 3 ) 

translation  to  go  from  a  SRT 
to  a  ORR  system 

TFM 

- 

k'T. 

:  % 

number 

Index  of  the  satellite 

CONTROL 

- 

- 

'  N 

number 

Index  of  the  time  point 

CONTROL 

No  ou  tpu  t  var  I  ab  le 


see  r  I  gure  1  2 
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gure  14:  Mechanical  coverage  area 
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4j  tp  I  *• 

Var 1  ab les 

"Mansions  | 

t  | 

Nature 

. .  ....  .  .  J 

Dest 1 nat 1  on  j 
prog' am 

Sys  tem  of  i 

i 

Coordinates  j 

Units 

n  jrnber 

! 

angle  from  the  grazing  angle 
tor  t  ie  Inside  border  o'  the 
coverage  area** 

STATE 

1 

radl  an 

GAN  +  ( G  ) 

number 

angle  from  the  grazing  angle 
for  the  outside  border  of  the 
coverage  area** 

STATE 

rad!  an 

XRC*(G) 

vector  (  3 ) 

distance  between  the  center  of j 
the  Earth  and  the  planes 
containing  the  boundaries  of 
the  coverage  area* 

TDE 

k  7i. 

RRC*(G) 

vector  (3) 

radius  of  the  boundaries  of 
the  coverage  area 

TDE 

- 

km. 

OT\<*(G)l 

number 

distance  between  the  satel  llte 

and  the  outside  border  of  the 

coverage  area 

STATE 

km. 

**  ,ee  r  I  gure  1  2 


see  r  I  gu re  14 


E 


r>~M 


"jocU.m:  calculates  soth  of  the  par  abaters  used  In  the  calculation  of 

the  position  of  the  borders  of  the  coverage  areas  for  the 
machanlcal  scanning  (see  Figures  12  and  14) 


1  nput 
Variables 

D 1  mens  1 ons 

Nature 

Source 

program 

System  of 
Coordl  nates 

Units 

•'v  ;(<;> 

O.NS/JT) 

position  of  the  satellite 

FOSSA T 

ERF 

km. 

V.'HO) 

number 

angle  from  the  gracing  angle 
for  the  outside  border  of  the 

coverage  area 

CONTROL 

~~ 

radl  ar 

R  (G ) 

number 

radtus  of  the  Earth 

CONTROL 

- 

km. 

T',  (o> 

number 

number  of  satellites 

- 

- 

■'ll  (G ) 

number 

number  of  time  points 

- 

- 

IS 

number 

Index  of  the  satellite 

CONTROL 

- 

- 

IN 

number 

Index  of  the  time  point 

CONTROL 

- 

- 

-  no  1 1  on: 

transforms  a 

3E0  system  Into  a  ESR  system. 

1  npu  t 
Var  i  ab  1 

0  1  mens  1  ons 

Nat-ure 
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program 

System  of 
Coord!  nates 

Uni  ts 

■*;(!/>> 

N*  x  1 

coordinate  of  a  point  in  the  CEO 
system 

BAR 

_ i 

ESR 

km: 

degrees 

ro'osPH 


Output 
Varl ables 


01  mans  Ions 


Nature 


Destl nation 
program 


System  of 
Coord I  nates 


Units 


MRC(LR) 


N*  x  3 


coordinate  of  a  point  In  the 
ESR  system 


km; 

radl ans 


**l  =  sny  Integer  number  greater  or  equal  to  1. 


jnetl  on: 


calculates  Hie  netrix  jf  rotation  needed  to  transform  a  cS’-i 
into  a  ERE  and  the  reverse. 


1  np  u  *■ 

/ar 1 ab les 

O I  nans  1  ons 

Nature 

I 

Source  i 

program 

i 

'  System  of 
Ooordi  nates 
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PO(LD) 

vector  ( 3 ) 
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POSSAT 

ESR 

km, 

radl ans 

VT(LD) 

number 

t  Ime 
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TOPO 


Output 
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MRO(LD) 


XO 

XOR.YDR 

ZDRtG) 

ZD 
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3x3 


Nature 


rotation  matrix  to  go  from  ESR 

to  ERE 


translation  to  go  from  ERT  to 
ORR 


Destl nation 
program 

POSSAT 

POSSAT 

POSSAT 


System  of 
Ooordl  nates 
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ZRf 

ERE 


POSSAT 


ERI 


PSX(LR)  N*  x  3  spherical  coordinates  PRSE  spherical  km.. 


P 


PRSM 


radians 


LvTT7. 


lornntAs 


Function:  transforms  spherical  coordinates  (in  an  array) 
into  rectangjlar  coordinates  (see  Fig.  2) 
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program 
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spherical  coordinates 
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PRSM 
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System  of 
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rectangular  coordinates 

BAR 

PRSE 
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km 

PRSM 

*N  is  an  integer  ranter  greater  or  equal  to  1. 
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rvKCsP 


met  1  on;  transf  or  in,  rectan  gul  ar  coordinates  (In  a  vector)  Into  spherical 


V  i'- I  ah les 

rsRti.D  > 

AREP(G) 


c oordl  nates 

(see  r  1  g.  2)  • 

01  mans  Ions 
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Sys  tem  of 

program 
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rvRCSP 
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jDestlna+Ion  System  of 
program  Coordinates 


l’S<(LR )  I  vector  (3)  vector  In  spherical  coordinates  BAR 


r  jncH  on: 


transforms  spherical  coordinates  (In  a  vector)  Intj  rectangular 
coord! nates  (see  -  I g.  2). 
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Nature 

Destl natl an 
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CONTROL 
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Coordl  nates 


system  of 
Coordl  nates 


jo :  ft  on. 
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to  the  SXT  system 
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PO(G) 

vector  ( 3 ) 

satel 1 Ite  position 

TWCSP 

ESF 

km. , 

r  adl  an 

MRS (G ) 

(IS,  IN,  3,3) 

matrix  of  rotation  to  go  from 

the  SXT  to  the  ERF  system 

POSSAT 

VOPOC 


lutput 

V.j r  i  ables 

R  1  liens  1  ons 

Nature 

Destf nati on 

program 

System  of 
CoordI  nates 

Units 

MPT (LR ) 

N  x  3 

coordinate  of  the  final  vector 

PRSE 

SXT 

km. 

j  i ;  r »  on: 


trans^orins  a  >^etor  MCR  from  the  system 
to  the  >XT  system 


N  -  n jnber  of  line*  of  f'ne  array  MC^ 
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System  o- 
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WSM 

::*r  i 

'n. 

VPS  (LD  ) 
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WSM 

i:r 

m. 

PO(G) 

vector  ( 3 ) 

sate  1 1 Ite  pos 1 tlon 

TVRCSP 
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k  m.  f 

radl  an? 

MRSE(G) 

( IS, IN, 3|3) 

matrix  of  rotation  to  go  from 
the  SXE  to  the  ERF  system 

POSSAf 

VOPOE 


")u  tpu  t 

Var  I  ab  1  ?s 

G|  -rens  Ions 

Mature 

Ges  1 1 nat I  on | 
pr  ograin 

System  o‘  j 
Go'Xd!  nati* 3 

‘ Ini  *s 

VPH'..R  ) 

M  x  3 

coordinate  of  the  final  vector 

PRSM 

G'T 

!  n. 

Conclusions 


We  described  two  complementary  software  packages  designed  to 
perform  studies  of  time  of  revisits  of  satellite  system.  Input  and  output 
variables  of  every  routine  were  described,  examples  of  the  output  were  given 
and  a  block  diagram  of  the  CONTROL  program  of  each  package  was  presented. 
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The  workspace  TSATCO  contains  the  following  programs:  The  DRAW 
package  of  APL  plus,  ANGR,  ANGROT,  ANGRV,  BROWNS,  CONTROL,  COURBE,  EARTH, 
HISTO,  0P0,  POSSAT,  PRSE,  PRSM,  RAX,  RAY,  RAZ,  STATE,  STATM,  SYSVEC,  TDE , 
TFF ,  TFM,  TGEOSPH,  TOPO,  TRECSPX,  TSPHREC ,  TVRCSP,  TVSPRC,  VOPOC,  VOPOE. 

The  workspace  SATCO  contains  the  same  programs  as  the  workspace 
TSATCO  except  that  the  program  CONTROL,  STATE  and  STATM  differs  from  the 
programs  having  the  same  name  in  SATCO. 


TSATCO 


VANGROTbJ]! 

V  VZ  ANGROT  l 7X-.AREP 
Cl]  +((*/(.VX*VZ))*d)/l  50 
[2  ]  *(  (  +/(VZ*VZ)  )*1  )/£51 

[31  -((+/( r/ATx  VY)  )*l  )/£,52 

[4  ]  PV**7Z[l];C'Vr*m2];C,.VZ*t'Z[3];C,4ir*mi];C'4i,*m2];P4Z*ia[  3  ] 

[5]  C3X<-(CVY*CAZ)  -(GAYxCNZ)  iCBY*-(CAX*CNZ)-(.CNX*CAZ)  ;CBZ+-(CNX* 
CAY) ~(CAX*CNY) 

Co]  VY+-CBX  ,CBY  ,C3Z 

[7]  *(( \CNZ)=l)/L2Q 

[8]  *(( |?4Z)=l)/£30 

C  9 ]  *( ( |CSZ)  =  l)/£40 

[10]  +(CBZ  =  0  )/no 

C  1 1  J  *(  £  4  Z  =  0  )  /  £3  0 

C  12  ]  VXP*-VXPi  (  (  +/  (VXPxVXP)  )*0  .  5)  -,V  XP*-V  SX*  {C  BY  ,C  BX  .0)  ;  (7SX*(  (  1  1 

0)x(/743r>O))+(C"l  1  0  )x  (CAY<0  )  ) 

C  13  ]  DEN0+-(  (CBX*C8X)+{CBY*C3Y)  )*0. 5 

C  1 4  ]  TYP+-AREP  \CBX  ANGR  CBY  iTXP-AREP  -.VXPLDANGR  VXPl  2] 

C  15  ]  SZ<-(TYP>TXP)*(  1-SZZ)  SZZ*  (  ?yp>  (  3  *  2  )  )  a  (  TYPS  (  2  *PI )  )*(TXP$( 

PH  2  )  ) 

Cis]  sz*  ( (  t^p>  ( 3  xpr  i-  2 ) )  a  ( rjf  ( 2  xp.r ) )  a  ( ryps  (  pi  *  2 ) ) )  +  sz 

C  17  ]  P«fy«-(  +  /(  (5Z=1)  ,-(SZ  =  0  )  )  )x  (  (S4Z>0  )-(S4Z<0  )  )*PHY \PHY+-~ 2o( +/ ( 

/^x/YP) ) 

C  18  ]  PH  X*-  (,  {  A  REP  x  (  C7BZ<0)-(.?BZ20  ) )  )x(SZ=l )  )  +  (  (  (?8Z<0  )-(?SZ>0  )  )x  ( 
Pr-4PPP)x(9Z=0)) iDENO  ANGRi \CBZ) 

C  19  ]  *£3  0  ;P7Z<--4ffSP;  ( tftfPC  1  ]  )  4V5P(  /JTPC  2  ]  ) 

[20  ]  £7  0  :DEN0*-{  UC  AX*C  AX)  +  (C  AY*C  AY)  )  *  0  .  S) 

[21  ]  PHY*-AREP*  (  xS4  Z )  x  (  (CBX'zQ  )  -  (CBX<0  )  )  x  (  (  S  A  If  £0  )  -  (  C  4  Y>  0  )  )  ;  DENO 
ANGR(  | C  4  Z ) 

[22  ]  PHX+PI*  (  (  (<7VZ>0  )x0  )+(SVZ<0  )  ) 

[23]  *£S0  ;P7Z<--47?SP;S4X  ANGR  CAY 

C  2 4  ]  £80  :  PHY*-0 

C  2  5  ]  DENO*-  (  (  (.CBXxCBX)  +  (CBY*C3Y)  )  *0. 5  )x  (  *  /  (  (SBZ>0  )  ,-(CBZ<  0  )  )  ) 

[25]  PHX*-- AREP’.DENO  A  NGR  CBZ 
[27]  *£50  ;P3Z^-4PPP;/74y  ,4<VSff  £4/ 

[28  1  £20  :  PHY*- (  0x(SVZ=l  )  )  +  (  9I*(CN Z-~l  )  )’,PHX*-  0 

[29]  *£50;P7Z<--(4/?5P-(PI»2));S9X  ANGR  CBY 

[30]  £30  :PHY*-(  (  P  T  *2  )  *  (C  AZ=  l )  )-(  (Pr*2)x(S4Z=  1)) 

[31]  ~*LoD  \?HX*-Q  \PHZ*-{PI  *2)  -  AREP  \  C8X  ANGR  CBY 

[32  ]  £4  0  :P.7Z*0  ;P7*<-(  (Pr*2)x(S9Z="l  )  )-(  (  PT*  2  )  x  ( S  BZ  =  1  )  ) 

[33]  *£30;P75f*4/?SPx(  {CBZ-  1)  ~{CBZ=D)  iCAX  ANGR  CAY 

[34]  £50 :  *  £9  1 ; *  ERROR  FROV  ANGROT  VZ  AND  VX  ARE  HOT  PERPENDICULAR 

f 

[35]  £5 1 : *£9 1 ; ' ERROR  PRON  ANGROT  VZ  IS  NOT  UNITARY' 

[35]  £5  2 : *£9 1 ; ' ERROR  FRON  ANGROT  VX  IS  NOT  UNITARY' 

[37]  £50  :  r7If*8'/,?;  xp^f  .  xPZ)  ;P41r  PHY\RAX  PHX\RAZ  PHZ 

[33]  *(  ( I /7[l ])50.9999)/ £90; VT-VX+ . xjyp 

[3  3]  *(  (  I  7  7C3  ]  )S0 . 99  9  9  )/£90  ;  VT<-VZ+.  *IVR 

[4  3  1  *£3  1  ;  *(  (  |  VTZ2  ]  )S0 . 99  9  9  )  /  £9  0  -,VT*-VY  +  .  *IHR 


TSATCO 


VA/YGROTCOl}^ 

[41]  £90  :-*-£92  ; '  RRROR  FROM  AN  GROT  TV  THR  ROTATION  ANGLRS 

[42]  L 92 : *  VX='  \VX 

[43 ]  » yy=*  xvy 

[44 ]  * VZ= • \VZ 

[45 ]  'PHY=' ; PHY 
[48]  'PRX='iPHX 
[47]  ,PVZ=»;P¥Z 


Cl  1 
[2] 
[3  ] 

[4] 

[5] 
[5] 
[7  ] 
[8  ] 
[9  ] 

C  19  ] 

cm 

C  12  ] 
C  13  ] 

cm 

C  15  ] 
C  1  5  ] 
C  17  ] 
C  1  8  ] 
C  19  ] 

C  2  0  ] 
C  2  1  ] 
C  2  2  ] 
C  23  J 
C  2  4  ] 
C  2  5  ] 
L25  ] 
C  2  7  ] 
C  2  8  ] 
C  2  9  ] 
C  30  ] 
C  3  1  ] 
C  3  2  ] 

C  3  3  ] 
C  3  4  ] 
C  3  5  ] 
C  36] 

C  37  ] 
C  38  ] 
C  39  J 

C49  ] 
C  4  1  ] 

C  42  ] 


TSATCO 

NCONTROLL  J]N 

CONTROL;  PERI ;  APOiLONR;  LA  TRi  9;  £;  RSS-,  XE\  YE-,  ZE'.CORi I;  /?;  Z£2  ; 

SC F ; TV Pi TOPiTOPiNLIMiNCAL; I HR ;IV 
' ST  ART  I N  G  TIME  OF  THE  SIMULATION-' 

TO+l 

'TIME  INCREMENT-' 

TI+3 

'DURATION  OF  SIMULATION  - ' 

TM<-T0+?D->-3 

'POSITION  OF  THE  OBSERV ATION  POINT' 

'ALTITUDE  ABOVE  EARTH  SU  RF  ACE-0 ' 

ROB-R+HR*- 0  -5  37  0 

’ LONGITUDE^ O' 

THOB+LON  R*PI  i  180  ;  LONR*-  0 
*  LATIT'JDE=0 » 

PHOB*-(.PIt  2  )  -  (  LATRxPIiiaO  )  ;  LATR+D 
PO-ROB ,TH0B,PH0B 
NS-lt (pSYSAT) 

OME-D  .  00437  iPREFl  ;3  ]•*-((  i  10  )-l  )*PH  18  ;PREF*-(  10  3)p6  370  0  0 

V957«-10p0  ;Zr2"-0  ;TREV*-ldp  (.TO- (  3xTI )  +  l )  ;  TREVS-(.  NS  ,  10  )  pO 

REVl4-REV2<-REV34-REVH*-REV5*-REV&*-REVl+REVB<-REV,i*-REVQ*-QpQ 

BOXxNINDOW  0  6.2  0  6.2  \STRAPIS  7  \SC  A  LS  (  (  -  SC  F)  .SC  F ,  (  -  SC  F)  , 

SCF)  \SCF*-{SY5ATl\  ;1  ]+SyS4T[l  ;2]  +  (  2*6370  )  )*2 

EARTH ;  MRF*-OPO  PO 

POR-TVSPRC  PO  i  MTO+-MRF 

NO  A  L+TD*TIiNLDt+9iTDP*-TDiINP+liT0P+T0iTVP+TMiT+T0iIN+l!tll<-l 
L  30  :-*(  NCAL>NLIM)  /  L20 
TMP+-TM ;  TOP*-T \  TDP*-NC  A  L*TI 

-*£2  5  ;  'TOP-'  ;  TOP  \  '  TMP=  '  \  TMP  iPOSSAT  (TOP  ,T  I  tTDP) 

£20  :  T0P*-T 

POSSATiTOP.TI  ,TDP)  ;  TMP*-TOP+(NLIM*  TI)  \  TDP-N  LIM*T  I 
£25  :  T*-TOP\  IN*-l  ;IS«-1 
£15  :-(rV/V=l)/£l7 
-*U/V=l)/£50 

*(PSO  Cl  ;  IS  i  TV]s(  -  P'2  Cl  3 ) )  /  £9  1 

*£17  ;0  DRA  W  { PSOL  3  ;  15  ;  IN-  1  ]  ,PS0[3;IS;  PV]  )  75 ( PSO C 2  ;  PS;  IV- 1  ]  . 

pso C2 ; rs; rv] ) 

£5  0  :-*(  PSO  Cl  ;  T9;  rv]s(  -P9[l]))/C91 

0  9ff4V(PS0ZCrS] ,PS0[3 ;  IS;  IN]  )V  S (PSO Yl I S ]  ,  PSO  [2  ;  T.S;  TV]  ) 

£17  :-*(SrS4Z,C  PS; 10 ]=0 ) / £50 

-*£91  \STATM\TFM\  ALD+SYSATl  IS;  11  ]  ;  VS54  N*-SYSATl  IS ;  1  2  ]  ;  PSS[  ;  £5 
IN]PRSM  PREF 

£6  0 : ST  AT Ei TF E i A LD+ SY S ATI  I S ; \ 1 ] ; PSSC ; 19; IN]PRSE  PREF 
£91  :-*(  (rS«-rS*l)SWS)/£15 

£11  :-*(  (T<-T+TI)ZTMP)  /£15  r='  ;  £+  n  ;  I/V-TV*  1  ;  IS«-1  ;  INN-INN*  1 

+(T>TM) / £35 
£VP*-£/VP*  1 

-*£30  ;.VC4£<-V.71£-(/V£r'/+l  )  ;  PSO  Z«-PS0  C  3  ;  ;  (  TV-1  )  ]  ;  PSO  Y<- PSO  C  2  ;  ;  ( 

r.v-i )  1 


TSATCO 


VCOURBEim 
7  COURBE  VV 

[1]  VW+iVViSTRAPIS  5 

[2]  L\+VVlVV\.lU\L2<rVViVWit>VW^ 

[3]  L1«-(L1*(L1<0))+(0xU1>0)) 

[4]  AXES-,BOX-,SCALE  1  10  ,  £1 .  L2  ;  WINDOW  0  S.  2  0  3  .2 

[5]  VX*-{  i(p/7))xl0*(p7lO 
C  s  3  o  DRAW  VV  VS  VX 

7 

VEARTRL^lV 

7  EARTH-,  ANG-,X-,Y\R-,HMAR‘,I 

[1  ]  ANG-*-{  0 , 1 100  )xoi* 50  -,HHAR4-R*20PREPl ;  3  ] ;  .<?<-3  37  0 

[2]  X-R*2oANG 

[3]  Y+RxloANG 

[4]  0  DRAW  Y  VS  X-,I*- 1 

[5]  0  DRAW {0 ,R)V  5(0 ,0) 

[6]  £50:-*-((I‘'-I  +  l)£10)/£5a;0  DRAW(  H  VARl  H  ,  H  V  ARl  I  ]  )  VSC  1 0  0 , 1 0  0  ) 

7 

7Hrsron37 

7  HISTO  VV 

[I]  V2-10p0 iSTRAPIS  5 

C  2  3  vw-ivv 

[33  /?5^-7/[7V[l]]-^7[7V[pl^i/]] 

[  '4  1  L2+-l\  +  INC-,  L  1-m  I/VC  1  3  3  ;  <7-*-l ;  INC-RGi  10 

[5  ]  £11  :♦(  (W+l)*10)/£ll;£2«-£2  +  r/ir£{£l«-£l  +  r.¥£;Y£[«n«-+/(  (77az;i )  a 

( 7/<£2  ) ) 

[6]  AXES-,  BOX-,  SCALE  0  10  0  ,  NMA  ;  NMA+NEl  p/YV.4  ]  ;  NHA-kNE ;  WINDOW  0 
6.2  0  6.2 

[7]  N*-\ 

[8]  £13  :-*■(  (AMf+l)S10)/£13;0  DRAW  VY  VS  VX-,VY*-0,NEiNl,NEi  /  J  ,  0  ; 
VX+-(N- 1  )  .  (  V-l  )  ,N  ,N 

[9]  • HISTOGRAM  VALUES=' ;NE 

[10]  'MEAN='  -,MEAN  ;VS4 //<-(  +  /NE)i(.pNE) 

[II]  'MEAN  SQUARE  ERROR='-,MSQ-,MSQ*-(*-/(.(.NE-MEAN)x(NE-MEAN)))i{? 
VE) 

7 

70P0[Q]7 

7  MRO-OPO  P0-,0ME-,THR-,VX-,VZ 

[1]  VX-VZ*-3pO 

[2]  7X[l]«-(  loP0[3])x2OP0[2];7y[2]<-(  1OP0[  3  ]  )  x  loP0[  2  ]  ;  VXi  3  ]  «-2  oPQ  [ 
3] 

[3]  mi>-(2oP0[3])*2OP0[2];irz[2]«-<2OP0[3])xlOP9[2];!rZ[3]«-13 
P0  [  3  ] 

[4]  VZ  ANGROT  VX 

[5]  MR0+-IHR 


TSATCO 


RPOSSATl J ] V 

7  POSSAT  PTiRi  LihVLiISiPSS 

[1  ]  HT+l  V5,VT)pVr;  VS<-14pSrS4r;/\rZ’*-p¥Z,;M7’^P7’[l  ]  ,  ( (  PTC  2  ]  *  i  L  ( PTC  3  ] 

pt[2 ] ) )  +  ?r[ i  ] ) 

[  2  ]  P*-S370  ;PSS«-PS0«-(  3  ,  VS  ,  VP)  pO 

L3  ]  4<-(  (2x/?)+srS47[  ;l]  +  SyS4'Z,[  ;  2  1 )  *  2 

[  4  ]  ff-  (  4  -  (  S  +  sy  S4  r  C  ;  1  ]  )  )  i  4  ;  MRSE*-  {VS.  NT,  3 , 3  )  p  0 

[5]  L-A*(1-{E*E)) 

[3  ]  F-BROUNS  VP;  I<-1  ;  MRO*-{  NT ,  3 , 3  )  pO  -.MTR*-{  XT  ,  NS  ,  3  )  p  0  ;  MRS*-  ( .VS  .  NT  , 
.  3)pO 

[7  ]  120  :-*■(  (  r«-Z>l  )<.NT)/L20  iMTRlh  ;  ]«-(  tfS , 3  )  p  (  XOR ,  YDR  ,  ZDR )  \t4R0lE  ; 

]«-Vr[l ; l^TOPO  PO 

[8  ]  RSS+-NL  *  (  1  +  (  NE*2oF )  )  ;.yL«-($(Vr,yS)p4)x(  1  -(VE*M£)  ) 

[ 9  J  rs<-l  ;PSP<-(.VZ\  3)  pO; 

[  10  ]  tio  :PS5[  ;2  ]•*■  ffSSCTSs  ]xloP[rS;  ];PS£[  ;1  ]«-ffSS[rS;  ]x2oF[  TS;  J 
[111  K/VZ  4/VG/?{m54(_44SyS4'Z’[rS;]));71VZ<-24(  7  +SYS/1'Z’[  TS  ;  ]  ) 

[12]  PSGC  ;TS;  ]«-$PS5i-.  *.¥/?;  r«-l 

[  13  ]  7,3  0  :  P50  [  ;  PS;  I  ]«-(  PSGl  ;  IS ;  H  -  NTRl  I ;  PS  ;])♦•.  x  NROl  I ;  ;  ] ;  MRSlISi  l 

;  1-NRiD’J+OPO  VPH+TVRCSP  PSGiiISiI'\ 

[  14  ]  VRSElISi  7;  ;  ]  <-MR\  VNZ  AN GROT { VGi  (  (  *  /  (  VS*  VG)  )  *0 . 5  )  )  ;  VS+-PSGI  ; 

rs ;  r  ] 

[15]  -(  (  r«-Z>l  )ZNT) /LOO 

[15]  -(  (  IS+IS+l  )S.V5) /L10 

7 


VPffSP[^]V 
V  7P5  PRSE  PR5 

[1]  PR  Si  ;2  3+-PPSC  ;  2  ]•*•(<?  Vffxj) 

[2]  PRS+TSPHREC  PRS 

[3]  P/?S«-7PS  /OP(77  PPS 

[4]  PSX*-TREGSPX  PRS 
[5  1  PflS-Pr-PS*[ ;3  ] 

7 


7PPSV[[]  ]7 
7  VPS  PRSM  PRS 

[1]  PffS[  ;2  ]-PP5[  ;2  l  +  COVffxD 

[2]  PRS*-TSPH  REG  PRS 

[3]  PRS+-V  PS  VO  PI  E  PRS 

[4]  PSX*-TRECSPX  PRS 

[5]  PDS+-PI  -  PSXl  ;  3  ] 


T5ATC0 


V/?4X[[]]7 
7  RAX  AR 

[1]  flJK-3  3 p  1  0  0  0.  (2o4ff)  ,  (-104/?)  .0,(104/?)  ,234/? 

V 

VRAYlW 
7  RAY  AR 

[1]  RY+  3  3  p  (  204/?)  ,  0,  ( 104/?)  .0  1  0  ,  (  -  104/?)  ,  0 , 2o4  R 

7 

7/?4Z[l]7 
7  RAZ  AR 

[1]  RZ- 3  3p(  204/?)  ,(-104/?)  ,0,(104/?)  ,(  204/?)  ,  0  0  0  1 
7 

7Sr42’/?[l]  ]7 

7  STMT/?;  VP4/V ;  7-94  V;  V RD DFO ;  NTP  i  VIRD\NJ\CA 
[1  1  VPAN*-10pPAfl  i  VGAR*-10p5AR  i  VRD*-ldpRDVAX 

[2]  0P<7-(PS*[  ;1  ]^//?D)  *(PDS<VGAV) 

[3]  NTP+-+  /  DFO 

[4]  -( YTP<1)/LS 

[5]  TDE 

[6]  VIHD-DFO/ (DFOxild)  iJ<-l 

[7]  LI  :  NJ-*-VINDlJ] 

C  a  3  CA+-'ZI2'l»FMT  VJ 

[  9  ]  -»{TRESZ(.S*TI)  )  /  Lo  -.TRES+-T- TREV  Si  IS ;  IM] 

[10]  TREVlVJh<-Ti  TRE<-T-TREVlH-J~\ 

[11  ]  NREVl  VJ] -/!//?£/[  7-7  ]  +  l  iTREV  SllS \RJ^T 

[12]  DUWY<-z'  REV'  ,  (74  [  1  ;  2  ]  ,  '  *-REV'  ,  CA  [  1  ;  2  ]  ,  »  ,  THE ' 

[13]  LS  :  +  ((J+-J+  1  )<.HTP) /LI 

[14]  L 5: 

7 

752’4T,V[G  ]7 

7  ST  ATM; VPAVi VG AN \VRD ; DFO ;NTP \V IN D ; NJ \C 4 
[1  ]  VPAH+-10pPAH  iVGAR*-l0pGA?1iVRD-13pRDVAX 

[  2  ]  DF0<-(PSXl  ;1  ]SVRD)  *(PUS<VGAM)  a(  (  (  (2*PI) -PSXL  ;  2  ]  )  S  (  VS?  A  H*  P I 
36 0  )  )  v (  PSXL  ;2’i*(WSCAN*PI*ZoO))) 

[3]  NTP+-+  !  DFO 

[4]  +(HTP<1) /L5 

[5]  rov 

[6]  VIRD+DFO/ {DFO*\  10)  ;,/«-l 

[7]  LI  :NJ<-VINDIJ] 

[3  ]  CA<-'  L12'  AFVT  A/J 

[9]  -*(TRESZ(  S*TI)  )/L5  -,TRES*-T- TREV  SLI S  \NJ  ] 

[10]  2,/?P7[Ven-Z,;'Z’/??<-2’-7/?fiT[  V«7] 

[  11  ]  V/?P^[Vt7]-*-V/?fi'7[  V-/]  +  l;fPS'75[f5;4/e7]^r 

[12  ]  DUWY+e'  REV'  ,CAll  ;2  ]  ,  ’«-/?5T'  ,C  4  [  1 ;  2  ]  ,  '  ,  THE ' 

[13]  LS  :  *(  )SNTP) /LI 

[ 14 ]  £5  : 


T  SAT  CO 


NSYSVECt<l~\  7 
7  SYSVEC 

[11  'INITIAL  VECTOR^' 

L  2  ]  /+-G 

[3]  'NUMBER  OF  VECTORS^'; 

[ '+  1  N+-rl 

[5]  'ANGLE  ( DEGREES )  BETWEEN  THE  VECTORS='; 

[51  4V-H 

[7]  'THE  ROTATION  15  AROUND  THE  AXIS  Z ' 

[8J  S[l;]«-t7;S<-(V,3)p0;I«-l 

[9]  LI  :-(  (I+-I  +  1  )<N)/L1  ;S[(  r  +  1  ) ;  3«-y+  .  xflz  ;  /?4Z  (  AN*  PI*  1*1 3  0  ) 

[10  1  S 
V 

NTDEir 1  ]  v 
7  TDE 

[  1  1  C0-*-2OxPrHdO;<7«-l  ;  r$,rS*-ZS2,S<-l  9p  0 

[2]  7,11  :COR<-(  XRClJI-HS)  .  RRCl  7]  .  0 

[3  1  rr«-i  jz^rsci  l-co/?^] ;  rsrsd >cos[2]  -,cor*-cor+  .*mroiin-,  ;  ] icon 

COR  +  TRSO',COR<-COR+  .  xVffSCTS;  TV;  ;  1 
[4  ]  L10  :ZTfl^flfl?[71xio4V2’;  YI’ff-'-flfltt  7]x2o4A?r;  ANT+CD*IT 

[5]  ,/1-tfS)  .YTR.ZTR 

[6  1  Z/?fS[r7’+l]«-i7<?/P[3l;y(;rS[rr4l]^0/?[2l;C0ff^^0/?+.xVff0[fV;  ;  1; 

COR-COR+TRSO-.COR+-COR+  .*MRSlIS-,  IN',  ;  1 
[7i  -(( rr«-rT*i  )si8)/tio 

[81  -*-£1 6  ;♦(  (  J+-J  *-i  )s2  )  /  £1 1  ;0  7?P4V  ZSTS  7S  rSTS 

[91  LIB :'ERROR  FROM  TFG  VX  AND  VI  4 RE  NOT  PERPENDICULAR' 

[ID]  7,15: 

7 

NTDML'}  ]  7 
7  TDM 

[11  "rff^ll  3p 0  ;</<-!  ;  PTR*-2  11  3pO 
[2  1  ANT-r{  (  i  11  )  -5  )*  WSC  AN*PIil  3  00 

[  3  ]  til  ;CT«[  ;2  ]<-, ,'?ffC*[  71  x  204  NT',  CTRL  ;  3  ]«-ffff?[, /lx  10.4  NT-,  CTRL  ;1  ]^K[ 

<71  -  HS 

['♦  1  PTRlJi  ;  ]<-(  (CTR+.  xVRSElIS',  7V;  ;  1  )+TRSO  )  +  .  x  y/?0[  TV;  ;  1 
[51  -(  (7<-7+l  )S2  )/511 

[  5  ]  Ztfr.S«-(  pm  1 ; ;  3  1 ) .  ( bPTRl  2  ;  ;  3  1 ) ,  (  ?m  1 ;  1 ;  3  1 ) 

[7  ]  y'5  7,S^-P2,ff[l  ;  ;  2  1  .  ( bPTRl  2  ;  ;  2  1  )  ,PTff[l  ;1  ;  2  1 

[ 8  J  0  DRAW  ZGTS  VS  YGTS 


TSATCO 


VTFE [□  ]  V 
7  TFE 

HS*-{+/(?SGl  ;  r  9;  rV]xPS5[  ;  rS;  rtf 1 )  )*0. 5 

PAN*-  10(ffx(io(5*Pr»S))tVS)  ;  7  4  V  <-  1 0  (  ffx  (  lo(  (  PT * 2  )  +  4£D )  )  f  H  3 )  ; 
ALD*-ALD*PH  130 
BRC*-ri5*3o(PAN  ,GAN) 

VRC*--30(PAN  ,GAN) 

XRC-<,  (  -BRCxVRC)*  (  (  (  ffxff)  +  ( VRC*VRC*R*R)- (  3/??x  Sff?) ) *0 . 5 ) ) *  C ( 1 
( VRC*VRC) ) ) 

RRC<-(VRCxXRC)  +  3RC‘,TRS0-PSGl  ;  r.5;  rtf  ]  -  V27?[  rtf ;  rS;  ] 
PPV4Jf^(((35-IPP[2])x(^5-^P^[2]))+(.‘?ffC[2]xpp^[2]))*0.5 
7 


V!TPV[i]]V 

7  tfv 

+/(PSS[  ;  r.5;rV]xpS5C  ;  I  7 ;  rtf]  )  )  *  0  .  5 
P4  V«-  lo(  Px  (  io(  5xpif5  )  )  f  35) ;S4.V-  lo(  R*  (  1  o(  (  PI  1 2  )  *  A  ZD )  )  J  75 )  ; 
4  LD*-  ALDxpi  i  1 3  0 
3/?'7«-'/Sx30(P4tf  ,;?4tf ) 

VRC*--3o(PAN  ,GAN) 

XRC*-{  (  -  BPCxMPO  +  C  (  (flxff)*(  '4RC*MRCxR*R)-{3RCx3RC)  )*0.  5)  )i  '  (  1 
(  VRC*VRC))  ) 

RRC*-(  VRC*XRC)  +  3RC’,TRS0*-UltZ)p(PSGL  ;  T  5 ;  I V  ]  -  MTRlI  N ;  T3;  ]  ) 
PDM4X^(  ((3S-XPC[2])x(PS-^?CC2]))  +  (ffPC,C2]xPff^[2]))*D.  5 
7 


7TGP0SP«n]7 
7  MRC*-TGEOSPH  VC 

VCL  ;2  ;2  ]xpii  180 

VffC[  ;3]<-(PIi2)-(  V^7[  ;3]xpr»183);  VRC*-MC 


7T0P3C037 

7  MRQ*-VT  T0P0  PO’.OVE 

THR*-P0l2l  +  (.0VE*VT)  \0VE*-  0 .  00437  5  ;  7X«-7Z«-3pO 
KO/MP^C  1  ]x(  loP3[3  ]  )*20THR) 

YDR*-(PQZl  ]x(  1OP0[3  ])xlor.¥ff) 

ZDP«-(P<?[1  ]*(  20P9C3  ]  )  )x  (  x7Z[3  ]  )  ;  7Z[  3  ]«-10P3[  3  ] 

7*[l  ]•*■(  10P'9[3])x20Z"/P;7Af[2]-*-(  1OP0[3  ]  )*\OTHR\  VXi  3  ]«-2OP0[  3  ] 
7Z[l ]♦-( 2OP0C3 ] ) x  2oTHR\  7Z[2]«~ ( 2oP£> [ 3 ]  )  x loTHR 
7Z  ANGROT  VX 
VR0*-IVR 


TSATCO 


7  IWimW  FSR-.MH 

[  1  ]  PSXl  ;  t  ]«-(  +/  (.FSRxPSR)  )  *0.  5  ;  PSX-(pFSR) p  0 
[2]  MH-( (FSflC j2]xpsff[ ;2])+(?Sff[;3]xF5ff[;3]))*0.5 
C 3  1  PS*C  ;2  3-‘-4/?PP-<-PSff[  ;2  3  4VS/?7  PS/?[;3] 

[4]  PSX[  ;3  3-4/?PP<-PS/?[  ;1  34VG/?r/  MR 

7 

7  mcc-tsphrec  mc 

ill  W7<?<-(  p  vs)  po 

[  2  3  MCC [  ;  1 1+-VCC  ;  1  1  x  ( loVC[  ;3  ]  )x  20MCI  ;2  ] 

C  3 ]  MCC C ; 2 VCC } 11* ( lO*?r ; 3 ] ) x io*f?r ; 2 ] 

[  '+]  W7S[  ;3  ]«-.Y£[  ;1  3x20,V<7[  ;3  ] 

7 


72,P/?<7 SP[  J  3 7 
7  PSX-TVRCSP  PS/? ;  MR 

[1  3  PSX[l]«-(  +/(PS/?xP5ff))*0. 5;PSY«-(pPSS)pO 

[2 ]  (FS/?[ 2  3xpSR[ 2  3 )+(PSff[l]x FSP[1 3 ) )*0. 5 

[3  3  PSXT  2  ]-*-4/?5P;  PSRl  1  ]  4VSP  PSPQ2  ] 

[4  3  PSJf[  3  3+-4/?PP;  PS/?[  3  3  4WS/?  MR 

7 


72’(?SPRS['337 
7  TRF—TV  SPRC  PQ 
[13  TRF-3  pO 

[  2  3  rffPC  l]-*-P0[l]x(  10PS[3  ]  )  x  20P9[2  ] 
[  3  3  7Pq,[2  3^PS[l3x(  loPS[3])xioP9C2] 

[ 4 1  TPP[33^PS[l]x2oPS[3] 

7 


VVOPOCl^V 

7  MPT -VPS  VO  PCS  MCR'.VXiVZi  PO 
[13  PO-TVRCSP  V?SiVX<-VZ-3p3iMPr-{pM?R)p') 

[2  1  'fP2’[  ;U^VfSp[  ;  1 3- trpsc  1  3 

[3  1  VPr[;2  ]<-*?/?[;  2  ]-7PS[2  3 

[4  1  '?PZ’[;3]^VSP[;3]-(/PS[33x(x7Z[3]));72[33^-loPS[33 

[5  1  MPT— MPT  *  .  x  (  flA//?S[  TS ;  TV ;  ;  3  ) 

7 


7f/SPSS[a37 

7  MPT-VPS  VQPOS  MCR\VX\VZ\PQ 
[13  PO-TVRCSP  VPS',VX-VZ-3pO  iMPT-(pMCR)  pO 
[2  3  VPT[;13^MCP[;13-1?PS[13 

[33  MPT[;23^VCP[;23-VPS[23 

[4  1  ¥PT’[;3  3^VC,1?[;3  ]-(7PS[33x(xr/Z[33));7Z[3  3-loPO[3] 

[5  1  MPT-MPT*  .*{$MRSRL  rS*,rV;;3) 

7 
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[1] 
[2] 
[3  ] 
[4] 

C  5  ] 
[6] 
[71 
[3] 

[9] 

[10] 
[111 
[12] 
[  13  ] 

[14] 
[15  ] 

[15] 
[17] 
[18  ] 
[19] 
[20  ] 
[21  ] 
[22  1 
[23  ] 
[24  j 
[25] 
[25] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34  ] 
[35  ] 
[  36  ] 
[37  ] 
[  38  ] 
[39] 

[>+0  ] 

[41] 

[42] 

[43] 

[44] 


SATCO 

VCONTROLZ 017 

COST  R0l\  PE  RI ;  APO ;  LON  R ;  LA  T  R ;  V;  S;  8;  £  ;  RSS ;  XE ;  YE-.  ZE ;  COR ;  I ;  R ; 

ZI 2  SC F  TMP  TO P  TD P N LI  M  NC  A  L  I N N  I N 
'STARTING  TIME  OF  THE  SIMULATION^' 

TO*  2 

’ TIME  INCREMENT  -  * 

rr«-] 

'DURATION  OF  SIMULATION =  * 

TM*TO*TD*l 

' POSITION  OF  THE  OBSERV  AT  ION  POINT ' 

’ ALTITUDE  ABOVE  EARTH  SURFACE-0 ' 

ROB*R*HR*-0  ;  /?  <-  3  3  7  0 
• LONOITUDE-O ' 

THOB*LONR*PI i  130  ;  LONR*0 
' LATITUDE-0 ' 

PH08*(PI  i2  )  -  (  LATR*PI *  130  )  ;  £4Z\-?«-0 

PO-ROB ,TH0B ,PH09 

NS-lt(pSYSAT) 

OME*0 . 00437  ;  P,??F[  ;3  >((  i  10  )-l  )xPm8  ’.PRSF*(  10  3)p6370  0  0 
/Vff^/^lOpO  ;zr 2*-0  -.TREVS*(NS  ,10  )  pTSS/ ;  TffS/^lOp  (70-  (5*11)41) 

RE  V 1  «-  RE  V  2  .?  /  3 -PS  / 4 -PS  P  5  -P£  /  5  - PS  /  7  -PS  / 8  -PS  7 9 -PS  / 0  -0  p  0 

MRF*OPO  PO 

POR+-TV SPRC  PO  ;  MTO*  MRF 

NCA  L*TD*TI’,  NLIM*  30  ;  TDP+-TD ;  I, VP-1  ;  TOP*TO ;  TMP*TM ;  I-  TO ;  I//-  7/VP-l 
£30  :*(NC4L>NLIM) /£29 
TMP*TM\  TOP*-T;  TDP+-NC ALxTI 

*L2S-,'TOP='  ;TOPi  '  7.VP=*  ;  MP;  POSSAT  (TOP.TI.TDP) 

L20:T0P*T 

POSSAT  (TOP  ,  TI ,  TDP )  ;  TMP*TOP*  (NLIM*TI)’,  TDP-NLIM*  TI 

£25 :T-T0P;rP-l ;IS-1 

L 15  :-*■(  INN=l)/Lll 

*(  r -V  =  1  )/£50 

>£17 

£,50  : 

£17  :-^(SyS4r[rS;10]  =  0)/£S0 

-►£91;  ST  ATM  \TFM\  ALD*SYSATlIS ;  1 1  ]  ;  VSC4P— SyS4r[  rS;12];PSS[  ;  15; 
I/V]PPSy  PREF 

£5  0  :  STATE’,  TFE\A  LD*SYSATiIS  ;  1 1  ]  ;  PSG[  ;  IS  ;  IV]  PPSS  PPSF 
£91  :-(  (  IS*- IS*  1  )SA£S)  /  £15 

Ll\i  +  ((T*T+TI)ZTMP)  /  LlS\IN*IN*l',IS*liINN*INN*l 
■*(  TlTM)  /  LOS 
INP*INP* 1 

-£3  0  iNC4L*NC4L-(NLIM+l)  -,PSOZ*PSOl3i  ;  (  I V  -  1  )  ]  ;  PS0Y*PS0  [  2  ;  ;  ( 

rv-i )] 

£35  :  *  PS/1='  ;PSF1 
'  REV2-'  ’.REV 2 
» PS/3=' ;PS73 
’  PSK4=’  ;PSK4 
* PS/5=’ ;PS75 


SATCO 


RGONTROL [Q45  ]7 
' REV 6  =  ' ;RSV S 

*  RSV1  -  ' 

’ /?,?y8=*  ;ff«:U8 

•  r?5/9=  *  ;flS79 
• REV  0  =  ’  ;7?5FO 
'Sr542’=,  •,5/547’ 
£ '  )VSTO  DAT  AGO 
£ ' ) S4/S  DAT  AGO 


B-13 


SATCO 
VST  4TB l ]]7 

7  ST  ATE; V  PAN ; VS  AN ;VRD ; DFO ; N  TP ;  VIND;NJ;CA 
[1  J  V  PAN*-\D  ^>P  AN  ;VS  AN-*-\O^S  AN ;  VR0<-\QpRDVAX 

[  2  ]  DFO  +  {PSXl  ;  1  ]<;  VRD )  a  (  PDS<  VS  A  V ) 

[3]  NTP*-*-  /  DFO 

[4]  -(VrP<l)/L5 

[5]  VIND*-DFO/ (DFO*i  10  )  ;J-l 

[5]  Cl  :NJ*-VINDL  <H 

[  7  ]  C4-*-'  ZI2  '  &FMT  NJ 

[  3  J  -+(TRESZ{5*TI))  /  LS  ;TRES*-T-TRSVS[15;NJ] 

[  9 ]  TREVl  NJ]*-T;TRE-T~TRSVlNJ] 

[10  ]  NRSVL  N-J  ]*-NREVL  N  J ] *  1 ; TREV SL  T  S ;  NJ  ~\*-T 
[  11  ]  DUMY*-<i'  REV'  ,CAl\  ;2  ]  .  '*-REV'  ,CAl\  ;2  ]  ,  '  ,TRE' 

[12]  C6  :-(  (J*-J+l  )$NTP) /LI 

[13]  •  ’ 

[14]  '  NTP-  '  ‘,NTP  ;  *  15='; IS;'  T='  ;T 

[15 ]  'DFO-' i DFO 

[IS]  ' TREVS=' ;TREVS 
[17]  ' TREV=' ;TREV 

[18  ]  '  N REV ~  '  ;NREV 

[19]  L 5: 

7 


7S74rV[ ]]7 

7  ST  ATM-,  VP  AN\  VS  AN;  VRD;  DFO  i  NTP ;  VIND  ;NJ;CA 
[1  ]  yP4.V«-10pP4V;t7<?4V-*-lOpS4.V;i7'?0«-lOp/?0'f4X 

[2  ]  DFO*-(?SXl  ;1  l^VRD)  *(.  PDS<VGAN )  *(  (  (  (  2xPJ)-PS*[  ;  2  ]  )  S  (  VS74/7*  PC 
36  0  )  )  v(PSJT[  ;  2  ]  5  (  VSS  AN  x  PI  f  3  5  0  )  )  ) 

[3]  NTP*-*-  /  DFO 

[4]  *(NTP< 1)/C5 

[5]  VIND+DFO/ (DFO*\  10);J*-1 
[5]  Cl  :.V-7<-7r/YD[7] 

[7]  CA**'  ZI7'  l\PUT  NJ 

[  8  ]  -(  TRESz  (  5  *  rD  )  /C5  ;2’P6'S^r-7,1'?5/S[  rs;/v-/] 

[9]  V./]^2’;7P5^2,-2,PS'7[.V/] 

[10  ]  VffmV7]^flpr/[V7]  +  l;2’PSTS[r5;/V7]«-C’ 

[  11  ]  DUWY*-e  'REV'  ,  0 A  [  1  ;  2  ]  ,  »  +-P97 «  ,C4[1  ; 2  1  ,TRE' 

[12]  C6:*(  (-W*l  )SVrP)/Cl 

[13]  •  ’ 

[14]  ' NTP- ' ; NTP ; '  TS=';IS ;•  T=';T 

[15 ]  » DFO-' ;DF0 

[16]  '  TREVS=' ;TREVS 

[17]  'TRSV=';TREV 

[18]  'NREV=';NREV 

[19]  C5 : 

7 


'D  3  F  1 

n  i  t  i  o  n 

i  F 

t  h  e 

v't  r 

i  h  l  '  ! 

h  V 

th  - 

1  ;ii>r 

)  ;  f  i  n  i  t  l  o  n  )  f  t  h  3 
r  •»  f  ?  r  i?  n  '  e  o  i  i  n  t  s 


TSATCO 

Flow  chart 

af  cgntrul 

(part  I ) 


>  ’  t  t  i  n  j  n  f 

3  i  r  3  n  3 1  ?  r  5 
in  packet 

1.33  tine 

involved 
firm  at  ion 

_ } 

r-  -  ■ 

L _ 

H  :  i  o  i  t 

0  .3  C  k  -3  t  ? 

the  l  a  a  t. 

1 1  :  ’  t  -ui  i  F  i  n  j  ex 

F) r  t  h  ;  p  a  c  k  e  t 
F  o  r  n  it  i  ) n 


3  3  t -up  of  index  to 
control  the  picket 
Fornnt  ion 


F  1 1  :  u  1  1 1  i  i  n  n  f  t  h 
i  i  t-  ->  l  l  l  t  c  position' 


Colculation  oF  the 
3  a  tell  it  3  positions 
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SATCO 


TSATOO 


VANGRllH 
7  AX  4, VS/?  \Y\VOA 
[1]  -(4jT*0)/£l 

[  2  ]  *12  \AREP*{  {PI  i2  )x  (  4  7>0  )  )  +  (  3  *PH2  )x  (  4  70  ) 

[  3  ]  £1  :V0A*{  (  4r>0  )  a(  AX>0  )  )  .  (  (  4  )  a  (  4  JlfO  )  )  ,  (  (  4  TO  )  a  (  1  *<  0  )  )  ,  (  ( 

4/<0)a(  4Jf>0)) 

C'O  AREP*AREPt  (PT- AREP)  ,  {Pl*AREP)  ,  (  (2x?r) -4/?£P)  ;  4PSP+-  3o|4y*4* 

[5]  AREP**/VCA*AREP 

[6]  £2: 

7 

74VSffmi7 

7  4  REP*  AX  ANGRV  AY  \V0  A‘,  LV  \V ARE 
Cl]  AREP*  ip  AX)  pO  ;  LV*pAX-tV*l  •,*{  {  p  AX  )  *  {  p  AY)  )  /  LI  d 
[  2  J  £11  :-(  AXl  V~\*0  )/£l 

[  3  ]  *L2\  ARE  Pi  VI*  (  {Pli2)*{  A  Jr[7]>0))  +  (3x  PI  r2)*(4Z[/JO) 

[4  ]  £1  :  VOA*{  (  4y[7]>0  )a(  U[  lf]>0  )  )  .  (  H/CHaO  )  a(  4^[/]<0  )  )  ,  (  (  AYlV}< 

0  )a(  AXl  7]<0  )  )  ,  (  (  4i,Cn<0  )a(  AXl  n>0  )  ) 

C  5  ]  4PS'<-4ff£.(Pr-4/?£)  APT*  ARE)  ,  (  (  2  *  PI )  -  A  RE  )  ;  4.?£<-~30|  4  Y  [  V  J  i  AXl  V  ] 

[5]  4/?£P[^W^'<?4x4/?5, 

[7]  £2  :*L12  ;-*■(  (  V*V*  1  )£  LV)  /  £11 

[8]  £10:’ ERROR  FROM  ANGRV  THE  ARGUMENT  DO  HOT  HAVE  THE  PROPER  £ 
ENGTH' 

[9]  £12: 

7 

VBRO<VN  SlQ  ]7 

7  F*BR0WN5  MT  i  NU G M‘,NB MRS ;  MB 

[1]  NU*G*M;G*2 . 4£  16;V—  5.97S2M- 

[2]  H8*(HUi A*A*A)  *0 . 5 

[31  MN9*(*{NT  tNS)  pVB 

[4  ]  MTP*MT-0.  01  ;  MTS*MT*0  . 01 lMT*MT* { ${N T ,NS) p5Y S ATI  ;9]  ) 

[  5  1  ME*${NT  ,R5)  pE’,  M9*MNB*MT \MBP*MTP*MNB \  M9S*MT5*  MR  3 

[5  ]  09*{  lo  MB)  i  {  (  l*{ME*ME)-{2*ME*2oMB)  )*0. 5 ) 

[7  ]  CBP*{  10MBP)»(  (  l*{ME*ME)-{2*ME*2oMBP)  )*0.5) 

[  8  ]  C9S*{  lo  MBS)*  (  (  l  +  (  ME*ME)  -(  2xtf£x  20  MBS)  )*0. 5  ) 
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